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I ' Imposing analytic properties to states and observables we construct a perturbative method to 

obtain a generalized biorthogonal system of eigenvalues and eigenvectors for quantum unstable 
y—i^ • systems. A decay process can be described using this generalized spectral decomposition, and the 

^ ' final generalized state is obtained. 

m '. 
■ 

O ; I. INTRODUCTION. 

Following previous ideas about the use of analytic continuation in the study of unstable quantum systems ||l| , , 
@5 0j ill papers Q and Q we have developed our own version, obtaining new results. But, as there were not a 
perturbative method adapted to our formalism, we were forced to study only solvable simple models. Therefore, this 

^ ^ paper is devoted to fill this gap introducing a perturbation method adapted to our formalism. Let us briefly review 

^ I ^ ] the problem and the history of the proposed solution. 

^ ■ The decay of unstable systems in quantum mechanics is usually described by a Hamiltonian of the form H = 
' Hfree + Hint, whcre Hfree has a discrete spectrum imbedded in the continuous spectrum, and H has a continuous 
^ ■ spectrum. In non equilibrium quantum statistical mechanics, the generator of the time evolution is the Liouvillc-Von 
^ Newmann operator which in some cases can be written as L = L free + ^int ■ There is an overlapping between the 
^ , continuous and the discrete part of the spectrum of L free , where in this case the discrete part of the spectrum contains 
the zero eigenvalue of the generator of the free evolution. It is expected that the infinite dimensional degeneration of 
^ . , the zero eigenvalue to be partially removed by the interaction, so that the zero eigenspace of L turns out to be the 
■ equilibrium subspace. 

E.g., for the Friedrichs model, a prototype of a decay system in quantum mechanics, with Hamiltonian 



>< 



H = n\l){l\+ duJLo\Lo){uj\+ dLoV^{\uo){l\ + \l){uj\} , G M+ , = , 
JO Jq 

it is possible to find conditions on the interaction for which the spectrum of H is M"*". The discrete part of 
the spectrum of Hq = ri|l)(l| + doj lu\oj) {uj\ is eliminated by the interaction. A complete set of generalized 
right (left) eigenvectors (('^-''''1) of the Hamiltonian H with eigenvalue uj € R"*" can be explicitly obtained, thus 
H = dco uj\u}^){uj^\. When Kj — > 0, tends to \uj) and therefore the discrete eigenvector |1) is not recovered. 

E.C.G.Sudarshan et al § obtained a generalized spectral decomposition for the Friedrichs model by considering a 
Hamiltonian obtained from the previous expression by replacing the integrals over by integrals over a curve F in 
the lower complex half plane. This extended Hamiltonian Hr is meaningful as the generator of the time evolution for 
state vectors $ for which (w|<i>) has a well defined analytic extension to the lower complex half plane. The spectrum 
of Hr is the curve F and a point zq located in the lower half plane between the curve F and the real axis. When 
T4; — > 0, zq tends to il, and the discrete part of the spectrum of the "free" Hamiltonian is recovered (see also reference 

§)• 
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The same complex spectral decomposition for the Friedrichs model was also obtained by T.Petrosky, I.Prigogine and 
S.Tasaki with a perturbative method using a "time ordering rule" by which an small imaginary part is included 
to avoid small denominators, with a different sign according to the type of transition involved. 

Latter, I.Antoniou et.al. ||9|], developed a formalism to deal with singular observables and generalized states, which 
they applied in the framework of the subdynamics theory. We gave our version of this formalism in papers |^ and . 
We also considered singular observables and generalized states to describe the time evolution of a quantum oscillator 
interacting with a field, describing the evolution of the system to equilibrium in the thermodynamic limit 

The analytic dilation method [|o| is also a well established procedure to determine the contributions of the 
singularities to the time evolution of unstable quantum systems. 

In this paper we consider states and observables with suitable analytic properties (as in papers and |^), such 
that the integrals over the continuous part of the spectrum can be deformed into integrals over curves in the complex 
plane. The physical motivation of these analytic properties is explained in In this way we can avoid the ill defined 
terms involving the differences between discrete and continuous eigenvalues of the unperturbed generator of the time 
evolution, and implement a well defined algorithm to obtain a complete biorthogonal set of generalized eigenvectors 
in one to one correspondence with the basis of the unperturbed problem. This biorthogonal system can be ued to 
give a description of the decay process. 

The paper is organized as follows: 

In section II we present the perturbative algorithm for the case of pure states, compute the complete set of right 
and left generalized eigenvectors and discuss the roll of the complex eigenvalues in the approximated expressions for 
the decay process. The method is applied to study the tunneling through a barrier in a one dimensional problem and 
to the Friedrichs model, were the results are compared with the exact solutions of reference In section III the case 
of generalized states and singular observables for the Friedrichs model is studied with the perturbative method. In the 
conclusions of section IV we summarize the results and discuss the physical interpretation of the complex generalized 
spectral decomposition. 



II. PURE STATES. 

Let us consider a system where the free Hamiltonian H has a continuous spectrum plus a single discrete eigenvalue, 

1 

and there is a general interaction H- 

1 

H = H + H 

r°° 

H =n\l)(l\+ dujLj\Lj){u;\ fieM+, 

^0 

1 ~" 



H= dLjV^\Lj)(l\+ dLjV^\l){uj\+ da; dLj'V^^,\uj}{u;'\, V^^,=V^,^, £ R+, (1) 



We suppose that the functions Vi^ and V^tj' are endowed with some analytic properties that we will specify below 

_ 
equation (|6|). The generahzed right (left) eigenvectors |1) and \uj) ((1| and (w|) of H form a complete biorthogonal 

set that we can use to expand any vector state |$), i.e. 



(1|1) = 1, {uj\u;') ^ 6{lj - uj'), = {uj\l} =0, 1$) = + / (2) 

Jo 

The amplitude of a state vector <I>, observed in the state vector 5', is given by 

/>OC 

(^1$) = + / duj{'i'\uj){uj\<i>), (3) 

Jo 

and therefore 



(*|$) = (*|/|$), duj\uj){uj\ 

Jo 





1 

If we try to construct a complete biorthogonal set of eigenvectors of the Hamiltonian H =H + H, depending 

1 

analytically on the small interaction parameter which we assume included in i7, the standard methods of perturbation 
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theory are not apphcable, because the superposition of the discrete and continuous part of the spectrum of H 
will produce non defined terms in the perturbation expansion. However, this problem can be avoided if we are 
satisfied to obtain information about terms of the form or exp(— ii/t)|$), where $ is a vector for 

which (p(ijj) — has a well defined analytic extension (p(z) to the lower half plane, and 4" is a vector for which 

tli{u>) = has a well defined analytic extension il^iz) to the upper half plane. In this case the functionals (z| and 

\z) are defined through the analytic expressions |4| ||^: 



(z|<f} - ^(z), {^z) = 
Then, the Cauchy theorem can be used in equation (^) to write 

(*|<I>) = + / = («'|1)(I|$) + / dz(*|z)(z|$), (4) 

^0 Jt 

where F is the curve in the lower complex half plane indicated in Fig. 1. Therefore, for the states <I> and observables 
^ endowed with the analyticity properties mentioned above, we can define the new "identity operator" 



Jr - j^dz\z){z\, m<^) = mm 



The functionals |1), \z) and {z\ satisfy 

(1|1) = 1, {z\z')^5r{z^z'), (I|z) = (z|I)=0, (5) 
where 5t{z — z') is the & distribution defined on the curve F (Jp dz'Sriz — z') f{z') = f{z)). 

1 

We can also define Hr =Hr + Hr, where 

Hr=n\l){l\ + f dzz\z){z\, Hr= f dzV,\z){l\ + /dzT?|I}(z|+ f dz f dz'V,,'\z){z'l (6) 
Jr Jr Jr Jr Jr 

where Vz, V- and V^z' are the analytic extensions to the lower half plane oi V^, V^j and V^^'. Therefore we postulate 

1 

that these functions can be analytically extended at least up to the curve F. The operators Hr and Hr verify 

11 

H 1$) = (^'1 Hr \^) and (^fl H |$) = (^'1 Hr \^) for the restricted class of vectors defined above, and therefore 
HrCSLH be considered as the generator of the time evolution. Being Q a real number, there is no superposition between 



the discrete and the continuous part of the spectrum of Hr ■ 

We can now proceed with a perturbation approach to the right eigenvalue problem 

Hr\<i>) ^ A|$). 

By assuming the expansion 

12 12 

I*) = I $) + I $) + I $) + , A=A + A + A + , 

1 

with respect to the small interaction parameter that we suppose contained in Hr, we obtain order by order the set of 
equations 

\ 

A-i/rj I $) =0 (7) 

\ 1 / 1 l\ 

A - i/rj I $) = {Hr - A j | $) (8) 

0\2 /l l\l 20 

A - i/r I $) = iJr - A I $)- A I $) (9) 



\ n / 1 l\ ri-1 2 n-2 n 

A - i/r I $) = Ltfr - A I $ )- A | $ ) A | $) (10) 
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Analogously, for the left eigenvalue problem we have 



and we obtain 



(^-li/r = A(*|, (*| = I + I + (I I + , A=A + A + A 



0/0 ON 

\ix-Hrj =0 (11) 

i/o o\ o/i l\ 

I ( A - i/rj = I iHr - \j (12) 

2/0 o\ 1/1 i\ 02 

I A - iJr = (* I ffr - A - (* I A (13) 



n-2 2 On 

I A - iJr = ( * I Lffr - A - ( * | A | A (14) 



It is convenient to define the projectors 

P<i = |l)(l|, Pr = j^dz\z){z\, Pz = \z){z\, (15) 

endowed with the following properties 

Pd + Pv=lT, PdPd = Pd, PrPr=Pr, P,P,, ^ 5t{z ~ z') P,, P^P^^P^P^ = Q 
PdHv =HrPd = ^Pd, PrHr=HrPr, Pz Hr^Hr Pz ^ z P,. (16) 
With all these equations in hand we can begin our computations. 

A. The discrete spectrum. 

Let us start computing the eigenvalues and eigenvectors of Hr corresponding to the unperturbed eigenvalue fi, so 
we take | <&) = |1} and A= f^- With this choice equation (|^) is verified. Then equation (^ gives 

J dz \z){n ~ z){z\ h) ^ I dz \z){z\ Hr A 

and therefore 



r Jr 



A= 0, Pr\ $) = Hr 

{n- Hr) 



Equation (S), gives 



and therefore 



r 



dz \z){fl - z)(z\ I) =Hr I $)- A 



211 2 1 11 

A- (1| Hr I $), Pr\ $) = —Pr Hr I 

{n- Hr) 



For I $} = |1) and A= ^1, equations (0)-(|T^) give no information about (1| $) (n — 1,2, ...). Thus, we impose the 

condition (1| $) = (n = 1, 2, ...), namely the usual choice in quantum perturbation theory. This choice will fix the 
normalization of the eigenvector. In summary, up to second order, we obtain the right eigenvector and eigenvalue 
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\^n) - |1) + I ^-^Pr Ht |1) + Hr | |1), 

Xn = n+{l\HT ^Pri/r|l). (17) 



Hr) 

Using the expressions for Hr and Hr given in equations ra), we obtain 



Xn = fl + / dz 



Vt- z 

If Vz and V- are assumed to be analytic functions of z in the lower half plane, the integral over F can be transformed 
into an integral over : 



r°° V V f°° ( /I 

Xn = n+ duj ^ " " =n+ dwl -inSiij -n)-v[ — - 
Jo il-uj + tO Jq I \uj-il 



= 0-^ diuv(^^j^j\V^\^ -nr\Vn\\ (18) 

where we see that the perturbation of the eigenvalue > gives, up to second order, a complex eigenvalue An with 
negative imaginary part (see Fig.l). Notice that the imaginary part would not appear if < 0. 

__. __. n 

Starting from (\I> | = (1| and A= ^, equations (p_l|)-(14j) and the imposed conditions (vj/ |1) = (n = 1, 2, ...) provide 

a well defined algorithm to obtain the left eigenvector. E.g. up to second order we obtain 

{^n\ = {l\ + {l\{Hr Pr ^-^] + (M i Hr Pr ^| , (19) 

{n-Hr)j V Hr) 

with the same eigenvalue An . As we are searching for a complete biorthogonal system, it is necessary at every step to 
define the normalized eigenvectors as 



l/o) - |$n)/v/(*o|*o), (/o| = (20) 
which satisfy the condition ~ 1- 



B. The continuous spectrum. 



Let us now compute the eigenvalues and eigenvectors of Hr corresponding to the continuous spectrum of Hr- So 
we now consider equations (|7|)-([To|) with | $) = \u) and A= m (u G F). Equation (|^) is satisfied, and equation (||) gives 



- z)(z| $) + (w - 17)|1)(1| $} = / dz\z){z\Hr \u) + Hr \u)- X / dz\z)5r{z - u), 

Jr Jr 

or equivalently 

{u- z){z\k) = {z\ Hr \u)-XSriz -u), {u - n){l\ = {l\ Hr \u). (21) 

1 

If we assume that Vzu is an analytic function of the two variables z and w, the term {z\ Hr \u) — Vzu does not include 
a term proportional to 5r{z — u), and therefore equations ( pT| ) give 

A=0 Pr\h) = ^ Q PrHr\u), Pd\h) = ^-^PdHr\u) (22) 

{u + i{]-Hr) {u- Hr) 

Equation (0) gives 
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r 2 2 f 11 112/" 

I dz \z)[u - z){z\ $) + (u - <$>)= J dz \z){z\ Hr\^) + |1)(1| Hr I A J dz \z)5r{z ~ u), 



and therefore 



(u-z)(z| I) = (2| Ar I A(5r(2-u), (w - I) = (1| | (23) 



Using (H) and (|2^) we can compute 



(^1 Hv\\) = I dz' 



r 



u + zO — z' M — r2 ' 

which does not include any factor proportional to Sy{z — u). Therefore, equations ( p3| ) give 

A=0 Pr|l> = ^ Q Pr iJr | 4), Pd| I) = ^Prf Ar | $) 

(u + iO-ffr) (u- Hr) 

In summary, up to second order, we obtain the right eigenvector 

I/O = 1^^) + I ^^^r Hr + ^-^Pd Hr | |^^) 

^{u + iO-Hr) {u- Hr) 

+ I ^-^Pr Hr + ^-^Pd Hr ) \u), (24) 

(u + iO-Hr) (u- Hr) 

with eigenvalue A„ = u. Notice that we have included the factor iO, namely we make the usual choice of outgoing 

solutions of scattering theory, to avoid the continuous-continuous resonances. 

, 



Equations (0)-(14) with | ^ {u\ and A= u G F are used to obtain the corresponding left eigenvector. Up to 



second order we obtain the left eigenvector 



{fu\ = {u\ + {u\ \ Hr Pr ^ „ + Hr Pd ^ 



' -1^1 " 

(u - iO- Hr) {u- Hr) 

{u\\HrPr ^-^+HrPd ^) ■ (25) 

(u-iO-Hr) (u-Hr), 



It can be easily verified that the vectors given in equations (^0|)-(^5|) satisfy, up to second order, the orthogonality 
relations 

{fu\fu') = Sriu-u'), {fn\fu)^{fu\fn)=0. (26) 
They also verify, up to second order, the "completeness relation" 



(vl/|$) = (vl/|/r|ci>), = J^du\u){u\^\fn){h\+ J^du\U){fu\. (27) 

For the operators H and Hr defined in equations (|^) and (^ we also have 

{^\H\^)^{^\Hr\<^), Hr^Xn\h){h\ + J^duX^\U){U\, (28) 

showing that the eigenvalues and eigenvectors obtained above provide the spectral decomposition of Hr- It should 
be stressed that the expressions for /p Sind H^ given in equations ( p7| ) and (^8|)are meaningful only when they are 
used to compute (*|Jr|$) or {^\Hr\^) for vectors ($) for which ^(w) = (wl^*) {(p{uj) = have a well defined 

analytic extension to the complex upper (lower) half plane. 
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C. The time evolution. 



Let us make a first application of the above results. The probability Pt of the state vector $t = e *^*<i>o at the 
time t of being measured in the state ^ is given by 

Pt = |(*|expHi/<)|$o)|' = |(«'|exp(-zi/r<)|$o)|', (29) 

where 

(4'|expHi/ri)|$o) =expHAoi)(*|/o)(/o|$o)+^d^^ cxpHut)(^'|/0(/„|$o). (30) 

If the exact expressions for the generalized eigenvalues and eigenvectors were known, equations (|2^) and (^ ) would 
give Pt exactly. An approximated expression will be obtained if eigenvalues and eigenvectors are given up to some 
finite order in the perturbation expansion. 

Let us consider the survival amplitude of the unstable state |1), given by 

{l\exp{-iHrt)\l) =exp(-zAoi)(l|/n)(/n|l) + J exp(-iut)(l|/„)(/,|l). 

1 

For a small coupling between the discrete and the continuous spectrum, i.e. if in the expression for H given in equation 
(|l|) there is a small multiplicative parameter e included in V^, and using the explicit expressions of the generalized 
eigenvectors given in sections IIIA and IIIB, we obtain 

(l|/n)(/n|l) = l + 0(£2), (l|/„)(/„|l) = 0(£2). (31) 

If Xn is computed up to second order from equation ([l8|), we have 

|(l|exp(-ii7i)|l)|2 = cxp(-27ryjt). (32) 

This expression is valid if ^ 1 and e'^t < 1, i.e. for small coupling between the discrete and the continuous 
spectrums and for not too large times. This is the standard approximated expression for a decay process which can be 
found in several text books on quantum mechanics |l^ . The deviations from exponential behavior (Zeno and Khalfin 
effects for shot and long times respectively), are not shown in this approximation. Concerning Zeno effect, notice that 



the time derivative for t = of expression (32), is of order e^. If we compute the survival probability keeping the 
second order terms coming from equations (|31[), we can obtain the Zeno effect, because the initial time derivative at 
t = turns out to be of order 

D. Tunneling through a barrier. 

We now consider some particular models. Let us first consider the one dimensional rectangular well problem with 
the following Hamiltonian, given in coordinate representation 

where 

V (x) = > " (34) 

. . o' . 

The spectral decomposition of H is well known. The discrete part of the spectrum has eigenvalues Ej be- 
tween and Vq. In coordinate representation, the corresponding normalized eigenvectors {x\j) are of the form 

exp ^Ty^^(K) ~ Ej)x^ , the minus (plus) sign corresponding to x > a {x < —a). Therefore the vector states {x\j) 

are concentrated inside the rectangular well. We are going to consider the simple case in which the parameters a 

and Vq verify the condition aW^H^ < where there is a single discrete eigenvalue Ei, and the corresponding 
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eigenvector is an even function of the coordinate x ||T^. There is also a continuous spectrum with generahzed 

eigenvalues Ek — j^k"^ + Vb (—00 < k < +00). The generalized eigenvectors \k) can be obtained as the solutions of 
the Lipmann-Schwinger equation. The vectors and |1) form an orthonormal basis 



\k') = S{k^k'), (1|1) = 1, (l|fc) = (fc|l) = 0, I =\1){1\ + J dk\k){k\. 



(35) 



We can now add to H the potential 



V (x) 



-Vi 



\x\ < b 
\x\ > 6, 



(36) 



with b ^ a and Q < Vi < Vq. Thus, a barrier of height Vq — Vi and length b — a appears. The potential V can be 
written in terms of the basis given in (p5|) 



where 



V^\l)Vii{l\ + / dk'\l)Vik'{k'\+ / dkVki\k){l\ + / dk / dk'\k)Vkk'{k'\, 



Vn 



dx{l\x){x\l) + 



+b 



dx{l\x){x\l) 



Vik = Vki 



-Vi 



dx{l\x){x\k) + 



+00 



dx{l\x){x\k) 



Vkk' = -Vi6{k - k') + Vkk' 



r+b 

Vkk' - Vi / dx{k\x){x\k'). 

J-b 



Due to the exponential decrease of {x\l) when x ±00, together with the assumption 6 ^ a, it can be expected that 
Vii and Vik be small. 

' 1' 01 
The total Hamiltonian H =H + V can be also given by H =H + H where 





H 



H = 



{Ei+Vn)\l){l\- 

+ OC 

dk'\l)Vik'{k'\ 



+ 00 



dk{Ek-Vi)\k){k\, 



dkVki\k){l\ 



dk 



dk'\k)Vkk'{k'\ 



Our interest is to describe the time evolution of the unstable state which is an even function in coordinate 
representation. Only the even part of the Hamiltonian is relevant to this process 





El even 



He 



{E^ + Vii)\l){l\ + 



dk{Ek-Vi)\ 

'^even ){k even \ ; 



+ 00 



dk'\l)Vik'{k',,,„\+ / dkVki\k 

even 



dk 



+00 



dk \keven)Vkk' {k^y^„ \ 



where \keven) = ^{\k) + \ — k)). If the integrals over are deformed into integrals over a curve F in the lower 

complex half plane, and if Vn and Vik are small, we can use the perturbation method of the beginning of this section. 



The correction to the discrete eigenvalue {Ei + Vu) of Heven due to the interaction is 



(1| Heven 



1 



Heven |1) 



+00 



El + Vn — Heven 

VikVki 



-iO 



{Ei + Vn)-i^ + V„-Vi 



-dk 



r' + oo 



VikVkiS 



{El + Vn) - + Vo - Vi) 



dk 



VikVkiV 



1 



{E, + Vn)-{^ + Vo-V,) 



dk. 
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If {El + Vii) > (Vb — Vi) the corrected eigenvalue has an imaginary part —i-^VilVli, where 



and the survival probability of the unstable state |1) is approximately given by 



27r/i 



E. Friedrichs model. 



We will finally check our method with the well known solvable Friedrichs model. This model can be obtained as a 

/n 1 

special case of equation p), with V^j^/ = 0. The Hamiltonian is H —H + H, where 

H=n\l){l\+ djjjuj\Lo){uo\, H= duj {V^\uj){l\ + V^\l){u;\} , 17 G R+ (37) 
Jo Jo 

As we considered in the previous section, the variable lo G can be transformed into a variable z G F, being F 
the curve in the lower complex half plane shown in Fig.l. The corresponding Hamiltonian iJp is given by 

Hr=Hr + Hr, Hr^n\l){l\+ j dzz\z){z\, Hr= J dzV,\z){l\ + J dz%\l){z\, (38) 

where Vz and Vy are the analytic extensions of and V^j . 

It is possible to obtain exact solutions for the right and left eigenvalue problems i?r|^) = and {^\Hr = A(\E'|, 
where |$) = + /p rfz z|z) and (^fj = + /p If for the sake of simplicity 

T](X) =X-n- / dz (39) 
Jr A - z 

verifies r]{Xn) = for just a single complex number Aji located in the region of the complex plane between and 
the curve F, the following spectral decomposition is obtained ^ 

Hr ^ \n\h){h\+ I duu\U){ful 



r 



where 



= / dz-^\z) 

V'/(Ao) I Jr Ao - z 

(7o| = ^i==((l|+ / dz-^{z\ 
a/??'(Ao) I Ao - z 

\U) ^\u)+ , + / dz- \z) 

' ' ' r,iu + ie) y ' ip (u + zO) -z ' ' 

(!u\ - («| + , ^" + / 7 \ (^1 ^ " e F. (40) 

These generalized eigenvectors form a complete biorthogonal system |l) , in the sense given in equations ( ^6| ) and (|2 
Using equations ([l^)-(^0|) we can obtain the approximated eigenvectors 
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with the eigenvalue 



V V— 1 

Xn = n+ I dz = n- dujV - \Vn\^ . (42) 

\l — Z Jr. LJ — iZ 



r 

Using equations (HJ) and ( p5|) we also obtain, up to second order 

l/u>-|«) + -^|l) + -^ / dz- \-Az) 

with eigenvalue = m (u e F). 

In order to compare these results with the approximate eigenvalues and eigenvectors obtained in equations (|4l|)-(^3|), 
we need an expansion of the exact expressions ( ^ ) in powers of the interaction parameter. We can use equation ( |39| ) 
to obtain the proper eigenvalue Aji, satisfying ??(Ao) = 0, from the iterative equation 

n+i f yJZ 

Xn=n + I dz , \ , ,71 = 0,1,2,..., x=n. 

-z 

For small values of the interaction parameter, the secuence Ao converge to a single solution An of ?y(Ao) = 0. Up to 
second order we have 

X,,^n+ f dz^^^. (44) 
Jr - z 

Up to second order in the interaction parameter, we also have 

-^^l-Udzj^^, ^^^ + 0{V% ^-^ = ^ + 0(^2) (45) 
\/v {Xn) 2 Jp [il-zj'^ Xn~z il-z r]{u±ie) u - il 

Replacing (^4|) and (|4^) in (^) we obtain the approximated expressions (^)-(^), so our approximated algorithm 
reproduces the results of the solvable model. 

III. GENERALIZED STATES AND SINGULAR OBSERVABLES. 

Up to now we have only considered pure states and observables which can be represented by projectors constructed 
with normalized vectors. The extension of the formalism to mixed states, represented by density operators of the 
form p = J2jPj\^j){^j\ ii^jl^j) = 1' Pj ^ 0' Sj Pj = ^^'^ observables of the same form (O = J^j is 
straightforward. However, singular observables as f^^ duj |Li;)(a;|, which are not operators of the form just mentioned, 
normally appear, and they can not be considered in the framework of the perturbation method presented in the 
previous section or its straightforward generalization. In this case, it is possible to consider the states as functionals 



□ O 0- These generahzed states 
. The corresponding perturbation 



acting on the set of observables represented by operators including singular terms 
are also useful to study the approach to equilibrium in the thermodynamic limit |^ 
method is therefore a non trivial generalization of the one in section II. This generalization will be the subject of this 
section and, as we will see, contains new and interesting features. 

Let us consider the Friedrichs model, presented in section II. The Hamiltonian is given by equation ( |37| ) and we 
assume, for simplicity, that ~V^. 

For this model we are going to consider the class of observables which can be represented by operators of the form 

/•OO /"OO /"OC /"OO /"OO 

= Oi|l)(l|+ / dLuO^\oj){uj\+ dio diu' O^^,\oj){uj'\+ duj 0^1 duj'Oi^>\l){uj'\ (46) 



where Oi — Oi, Oui — O^j, Ouiui' — Ow'w, and Ouii — Oiw, being Oui, Ouu' , Oui and Oiw' regular functions of the 
variables w and w'. The second term of the r.h.s. of the last equation is a singular operator. 
In the Heisemberg representation, the time evolution of an observable is given by ^ 



^Superoperators, i.e. operators acting on operators, will be written in blackboard bold types. 
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-ijOt=l.Ou LO=[i/,0], 

where L is the Liouville-Von Newmann superoperator. 

It is convenient to introduce the notation \0) = O and to define the generaUzed operators 

|1) = |1)(1|, H = |c.)(c.|, \tou^')^\u){Jl M) = M(1|, = (47) 

Therefore we can write 

pOO pOO pOO pOG pOG 

\0) = Oi\l)+ dujO^\uj)+ duj duj' O^^, \ujuj') + duj0^i\Ljl)+ duj'Oi^,\W). (48) 
Jo Jo Jo Jo Jo 

When the states of the system are represented by the usual density operators p = , the mean value of an observable 
O is given by {0)p = Tr{p^O), and any density operator has the following well known properties 



Tr(pt) = Tr{p^I) = 1, Trip^iaA + fiB] = aTr{p A) + l3Tr{p^B), Tr{p^ B^) = Tr{p^B), (49) 

where / is the identity operator (/ = |1}(1| + duj\uj){Lu\ ~ |1) + duj\uj)). 

As we are going to consider the possibility to have more general states than the usual density operators p, we will 
consider generalized states as functionals p acting on the space of observables defined in equation ( H) . The action of 
a state functional p on an observable O will be denoted by {p\0) and it gives the mean value (0)p of the observable 
in the state p {{0)p ^ {p\0)). For any state functional p we assume the following properties 



(p|/) = l, {p\aA + pB)^a{p\A)+p{p\B), {p\B^ ) = {p\B) , (50) 

which are the generalization of the properties for the density operators p given in equations ( ^9| ) . 

For the class of observables given in equation ( |4^ ) it is convenient to define the hmctionals (1|, ((jj|, (woj'I, {uj1\ and 
(Iw'l, in such a way that 

(l|0) = Oi, {uj\0)^0^, {ujuj'\0)^0^^., (wl|0) = 0^i, (lu;'|0) = Ow. (51) 
Using (ll), (|0|) and (|l|) we obtain 



(p|0) = (p|l)(l|0) + / du:{p\io){oj\0) ^ duj duj'{p\u;uj'){ujuj'\0) + 



OO 



daj{p\ujl){u;l\0) + / doj' {p\luj'){luj'\0). (52) 

^0 



From this equation we obtain 

{p\0) = (p|I|0), 

where I is the identity superoperator 



I=|l)(l|+ / daj\u;){uj\ + / duj duj'\u;tj'){ujuj'\+ / dLj\ujl){u;l\ + / daj'\lu;'){luj'\. (53) 
Jo Jo Jo Jo Jo 

The generalized state functionals and observables defined in (^^ and (^) satisfy the orthogonality conditions 

(1|1) = 1, iu;\io') = Siu;-u;'), 

(lie.) = im') = im) = im') - O, (^c^'ll) = {u;u;'\0 = {u;u;'m ^ {u;u;'\ia = 0, 

(c.111) = {ujiio = i^m') = (^i|in - 0, (W|i) = (Wio - {wm = a^'m = o. (54) 
1 

Defining the superoperators L and L by 
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L O = [H,0], 



hO= [H-.O], 



1 

=L + L, 



and using the generalized operators and functionals defined in equations and (|5l|), we have 



dLu'{n~Uj')\lu}'){lu}'\+ / duj{uj - n)\u}l){ujl\ + duj duj'ioj ~UJ')\UJLU'){UJUJ'\, 



K.|l)+ / di0'V^,\Lu'Lu) 



L 



L - / dcuV^ [\uj1) - \luj)] (1| + / dcuV^ ~ \oj1)] + / duj 



duj 



V^\l)- / dij'V^,\ujuj') 



{uj1\ + / diU / dij'[V^\luj')~V^,\LUl)]{LUCj'\. 



(55) 



(56) 



From equation (^5|) and the orthogonahty conditions (|5^), the generahzed states and observables defined in equations 
( pT| ) and (|4^) form a complete biorthogonal set of generalized eigenvectors of the " free" superoperator L 



L|l)=L|w)=0, l.\luj') = {n-Uj')\lu}'), h\ujl) ^ {UJ 1.\U}LJ')^{U}-U}')\UJUJ'), 

(1| L = (t^l L= 0, iW\l= {n~uj'){luj'\, {lu1\1^ {uj ~n){Lul\, {iOLu'\l= [uj -lu'){luuj'\, 



(57) 



If we try to construct a complete biorthogonal set of eigenvectors of the Liouville-Von Newmann superoperator 

1 1 

L =L + L, depending analytically on the small interaction parameter which we assume included in L, the standard 
methods of perturbation theory are not applicable: the superposition of the discrete and continuous parts of the 

spectrum of L will produce non defined terms in the perturbation expansion. 

However, this problem can be avoided if we are satisfied, as in the pure state case, to obtain information about terms 
of the form {p\0) or (/9|L|0), considering states p and observables O for which (pjcjcj'), {p\u!l), {p\luj'), {ujuj'\0), {oj1\0) 
and (Iw'jO) have well defined analytic extensions {p\zz'), (p\zl), {p\lz'), {zz'\0), {zl\0) and (lz'|0) to complex values 
of z {z') on the upper (lower) complex half plane. We also consider states p and observables O for which {p\uj) and 
(wjO) have well defined analytic extensions {p\z) and {z\0) for complex values of z near in the upper and in the 
lower complex half planes. The physical bases of these analyticity properties are given in papers Q and |^]. 

For these states and observables, the Cauchy theorem can be used in equation (B2|) to write 



{p\0) = (p|l)(l|0) + / dcj{p\Lu)iuj\0) + / dz _dz'ip\zz'){zz'\0) + 



dz{p\zl){zl\0) + _dz'{p\lz'){lz'\0), 
r Jr 



(58) 



where F is a curve in the lower complex half plane, as shown in Fig. 1, and F is the conjugate curve in the upper 
complex half plane. Therefore 



ip\0) = {p\le.t\0), W = |1)(1| + / da;\Lu)iuj\ + dz dz'\zz'){zz'\ + / dz\zl){zl\ + / dz'\lz'){lz'\. (59) 



If in the expression for text given in equation (p9|), the term ^ duj\uj){uj\ is replaced by ^Ydz\z)[z\ or by Jp (iz'|z')(z'|, 
the identity {p\0) = (pjlgictlO) is also verified for the restricted classes of states and observables with the well defined 
analytic extensions defined above. 

1 

We can also define l^ext =hext + hext, where 




1 



dz'{n- z')\lz'){lz'\+ / dz{z^n)\zl){zl\+ dz dz'{z- z')\zz'){zz'\, 



= / dzV,\zl){l\ - / dz'K.|lz')(l|+ / dz'V,>\lz'){z'\~ / dzV,\zl){z\ + 



dz'V^'\l){lz'\+ dz dz'V^\zz'){lz'\+ / dzV^\l){zl\ 



dz / dz'V^>\zz'){zl\+ _dz / dz'V^\lz'){zz'\- _dz / dz'V^'\zl){zz'\ 
r Jr Jr Jr 



(60) 



(61) 
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1 11 

The superoperators hext and hext verify {p\ L |0) = (p| hext \0) and {p\ L |0) = (p| Le^i |0). The generahzed 

eigenvalues of 'Lext are shown in Fig. 2. The form given in equations ( |60|) and (|6ll ) for the superoperator h^xt —Liext 
1 

+ 'hext can also be obtained directly from the definition h^xtO = 11^0 — OHr, where Hr was defined in section II. 
For the right eigenvalue problem Lea;(|$) — A|$) we assume an expansion 

12 12 

1$) = I $) + I $) + I $) + , A=A + A + A + , 

.... 1 

with respect to the small interaction parameter that we suppose it is contained in L- The following set of equations 
is obtained order by order 

\ 

X-hextj I $) = (62) 
\ 1 / 1 l\ 

X-hextjl^) ^ [hext- X]\'S>) (63) 
00\2 / 1 l\l 20 

A - hext] I $ ) = {hext - A j I $ )- A I $ ) (64) 
\ n /l l\ n-1 2 n-2 n 

A - hext I $ ) = hext - A I $ )- A I $ ) A I $ ) (65) 



For the left eigenvalue problem (^'ILg^jf = A(^|, the expansion 

12 12 

(^-l = (^f 1 + (5- 1 + (4- I + , A=A + A + A + 



2:ives 



/o \ 

( * I ( A - Lext j = (66) 
i/o o\ o/i l\ 

( * I ( A - Lext 1 = ( * I ( Le:rt - A 1 (67) 

2/0 0\ 1/1 l\ 02 

( * I ( A - Lext j = ( * I (Lext - A j - ( * I A (68) 

/O \ n-1 / 1 l\ n-2 2 On 

( * I A - Lext = ( * I hext - A - ( * |A (*|A (69) 



It is convenient to define the following projectors 

Po = |l)(l|+ / dLj\Lj){uj\, Ppp = [ dz [ dz'\zz'){zz'\, Pp^ / dz\zl){zl\, Pr = / dz'\lz'){lz'\. (70) 



Jv Jr Jr Jr 







These projectors commute with Lext and also satisfy 



'0 + ifrr +Jr'p +rr, 




^0 Lexf=Lext Po 0. 





Therefore Pq is the projector on the invariant space of the time evolution generated by the superoperator Lext- As in 
the pure state case we will compute the discrete and the continuous spectrum. 

A. The discrete spectrum. 

We first consider equations (pa)-(p3) with 
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A= 0, 



$) = Pol $) = |1)(1| $) + / duj\uj){uj\ $) 



Multiplying equation ( |63D by Pq and by Qo = iext — Po, we obtain 



(-1) 



10 

A I =Po hext Pol 



h\ = 



1 
'0 Lext 1 $)• 



(71) 



(72) 
(73) 



iO ibext VQ 



Equation (p3) gives 



2 

A I $) 



h\ 1) 



. ; (-1) 

- ILext 7, 



HO ll^ext iro 



h l.ext Qo 







(-1) 



^0 ILext VO 



1 1 



1$). 



(74) 
(75) 



From the definitions (|7^) and (|6l| ) we have Pq Lea:* Po = and therefore A= 0: the degeneration of the space expanded 
by the projector Po is not eliminated by first order corrections. The previous equations give no information about 

n n n 

Pol $) ("- = Ij 2, ...), and we make the usual choice | $) = Qo| $) {n ^ 1, 2, ...). 

But the degeneration is partially removed with the second order correction. Equation (fz^), an eigenvalue problem 

2 

for A, gives 



2t:iV^ |1)(1| $) - 2TriV^ \1){lu = n\ $) =A |l)(ll $)+ A J dio\Lo){uj\ $). 
This equation has the solutions 

Xd^2TriV^, |-I'd)-|l), L=0, \L)=S{Lj-n)\l) + \uj). 

We give in appendix A the details of these calculations. 
For the left eigenvalue problem with 





A= 0, 



(^| = (^ |Po-(*|i)(i| 



(itj(vl/ |ci;)(a;|, 



equations (66)-(|69[) give 



(*| A =(*|(PoLex*Po), 



1 (-1) 

= (* I L.ext Qo TT"^ 



2 



(* I A = (* I 



(* IQo = (* I 



Po Le 



Jo Le:Et Qo 
(-1) 



h he 



1 1 
Lext — A 



(-1) 



h hext Qo 



(76) 
(77) 

(78) 

(79) 
(80) 

(81) 
(82) 



1 

Once again, the degeneration is not removed by the first order corrections, and A= as a consequence of Pq hext 

n n n 

Pq = 0. As there is no information about (v|/ |Po for n > 1, the usual condition (v|/ | = (il/ |Qo can be imposed. For 
the eigenvalue problem (|8^) we obtain 

2TriV^ {I - 2mV,Uh'^)i^ = ^\ =X I f dij{l, \u;){lj\. (83) 
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This equation has the solutions 

L= 2niVi, {^d I = (1| - (cc; = r!|, L= 0, (L | = (84) 

(See appendix A for the details). 

The right and left eigenvectors given in equations (|7^) and ( p^ form a complete biorthogonal system for the Pq 
subspace, i.e. 

00 00 , 00 00 

(^-d I $d) = 1, (^-^ I $^') = S{uj - oj'), (fd I $^0 = (*^ I ^d) = 0, 

/■ 

Fo = I $d)(*d \+ duj\ $^)(^'^ |. 



The degeneration of the zero eigenvalue of Lea* have been partially removed by the interaction: a proper eigenvalue 
Xd = 27r«VQ with a single eigenvector is obtained. However, the infinite degeneration of the eigenvalue Xu, = 
{u! G M+) remains. This is a new feature that must be taken in consideration: the infinitely degenerated eigenvalue 
corresponds to the invariant states. As we see from equations (^4|), these invariant states are expanded, up to zero 



order, by the singular states (vfi^ | = {uj\. These invariant states only appear if we introduce the singular structure 
for states and observables. The singular structure is also necessary for the description of the decay of the unstable 
discrete state. As we will see below, the imaginary eigenvalue is associated with the exponential approximation of 
the decay process. 

B. The continuous spectrum. 



Let us consider the following generalized eigenvalues and eigenvectors of Lest 



Xui=u-V,, I $„i) = |iil), I = u e r, (85) 



Equation (B3) gives 



{u-n-lext)\h)^iLxt~x)\ui). (86) 



When this equation is multiplied from the left by {vl\ we obtain 

1 1 
{u - u)(wl| $) = - A (5p(u - v), 

where (5p(u — v) is the " (5-function" on the curve T (Jp(iw(5p(M — v) f{u) = /(«)). From this last equation we obtain 



Xiii = 0. The first order correction for the right eigenvector can be obtained from equation (86) 



I ki) = -^|l) - / du'4^\uu'). (87) 
From equation (^^ we obtain the first order correction of the left eigenvector with the same eigenvalue 

(*«i I = ^[(1| -("!]-/ du'4^{uu'\. (88) 
The second order correction to the eigenvalue can be obtained multiplying equation ( |64| ) by (i;l| from the left 

2 112 

[u - v){vl\ $„i) = {vl\ 'Lext I $)- A«i %(u - w) = 



2 



du^^- Xui 

p u' — il 



(5p(u-w + -. 

u — il 



From this equation we obtain 
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Y2 roc, y2 /-oo / 



r 



, _ 

Therefore, the generahzed eigenvalues of Leat, represented by points z = u — Vl (u (zT) oi the complex plane, are shifted 

2 

to the upper half plane by the interaction. The corresponding eigenvalues of J^ext are the points z = u — Q.+ Aui, 
where u & V and Aui is given by equation ( |89| ) (see Fig. 2). 



We can also consider the following generalized eigenvalues and eigenvectors of Lea;* 

\iu'^n-u', I = (!-!„, I = (im'I, w'er, (90) 



The perturbation method gives 



Mw =0, 

= (91) 

This result shows that the generalized eigenvalues of equation (|90| ) are also shifted to the upper half plane by the 
interaction (see Fig. 2). 



Finally, we can consider the following generalized eigenvalues and eigenvectors of Lext 

An«'= u-u', I = \uu'), {^uu' I = {uu'\, u e r. It' G F, (92) 

(see the shaded region in Fig.l). In this case the perturbed expansion gives 

1 2 
\uu' \uu' — 

|$w) = -\|W) + ^^M), 

u — \ l u — \l 

{huu'\ = ^{lu'\ + ^^{ul\. (93) 

There is no shift of this part of the generalized spectrum by the interaction. 

C. Time evolution. 

In the previous subsection we obtained by a perturbation method a set of generalized eigenvalues and eigenvectors 



of the superoperator 'hext , which reduce to the complete biorthogonal set of hext when the parameter of the interaction 
goes to zero. The generalized eigenvalues of hext are shown in Fig. 3. It seems reasonable to assume that the obtained 
generalized eigenvalues and eigenvectors of hext, also form a complete biorthogonal set, at least for small values of 
the interaction parameter If this is the case, the time evolution operator is 

exp{ihextt) = J exp(iA^t) \^^){-^^\ + exp{i\dt)\^d){'i'd\ + J_du exp(iA„ii)|«'«i)(*«i | + 

+ J du' exp{iXiu't)\^iu'){'ifiu'\+ J_du J du' exp(iAwi)|«'w )(*««'!. (94) 



^This fact can be explicitly verified up to second order for the generalized eigenvectors obtained in sections IIIA and IIIB. 
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Up to second order, the eigenvalues are 

A,. = 0, 



\xu' r2 - u' + ii^Vl - dw Vl V 



1 



1 



Ku' ^u-u' (95) 

and the approximated expressions for the eigenvectors are given in the previous subsection. 
The time evolution of the mean value of an observable O, given by 

{0)p, = {pt\0) = (po|cxp(^L^)|0), 

can be easily computed using (p^). Let us consider the time evolution of the unstable state given by (po| = 

If we call by e the small multiplicative interaction parameter included in V^, we can use the explicit approximated 

expressions of the eigenvectors to prove that 

1) -0(£2), 
l)-l + 0(£2), 

and therefore the survival probability of the unstable state which is represented by the state functional (pol = (1| 
(and also by the vector state |1}) is approximately given by 

{pt\l) = (1| exp(iLt)ll) ^ exp{-2TTV^t). (96) 

This expression for the survival probability is valid if ^ 1 and t < e^^. 
Taking into account that 



(i|$^0(*-'M 


= S{uj' - 




(i|$d)(*dM 


^ -6{ij - 




(il$,i)(«',iM 


= 0{e'), 




(i|$i„0(^'i«'l^) 


= Ois% 




(l|$w)(*wk) 


= 0{e'), 




= (1 we obtain 






) = (l|exp(iLt)|a); 


^ S{uj - 


n) [1- 



for the initial unstable state {po\ 

{pt\uj) = (1| exp(iLt)|a)) S{uj exp(-27rVf^i)J , (97) 

also valid if <C 1 and t < e^^. Equations (@)-(0) describes the transition of the unstable state (1|, into {uj = 

Therefore, the obtained set of generalized eigenvectors is a useful tool for a complete description of the decay 
process. The pure state |1) is given in this formalism by the state functional (1|, and the standard approximated 
expression its the survival probability is reobtained in equation (p6|). But in addition we obtain in equation (|9^ ) an 
explicit expression for the by-products of the decay process. 

IV. CONCLUSIONS. 

When faced with the problem of the spectral decomposition of the Hamiltonian 

1 

H = H + H, 

H du;u;\uj}{uj\, 17 £ M+ 

Jo 

^ pOC pOC pOO pOO 

H = dujV^\oj){l\ + / dujV[^\l){Lj\+ du diu'V^Mi^'l 
Jo Jo Jo Jo 
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we find that the superposition of the continuous and the discrete part of the spectrum of H makes impossible to apply 
the standard methods of perturbation theory. 

If we restrict the class of states (observables) of the system to the set represented by vectors $ (^P) for which (cj|<i>) 
((wj^)) has a well defined analytic extension {{z\^)) to the lower (upper) complex half plane, it is possible to 

decompose the amplitude of the state |$) to be observed in the state I^E") in the form 

(*|$) = (*|/r|$), /r = |l)(l|+^dz|z)(z|, (99) 

being F a curve in the lower complex half plane as shown in fig.l. 

For the amplitude of the time evolved state |$t) = exp{—iHt) |$) to be observed in |\1/), we also have 

{^\^t) = (^'1 expi-iHt) 1$) = (*| eM-iHrt) |$), (100) 

where 

1 

Ht = n\l){l\+ j dzz\z){z\, 

Ht - j dzV,\z){l\+ j dzV;\l){z\+ J dzj dz'V,,j\z){z'\. (101) 

Starting from the right (left) eigenvectors |1) and \z) ((1| and (z|), which form a complete biorthogonal set for 
. . . . . 

the spectral decomposition of Hv, it is possible to obtain a well defined perturbative algorithm to compute the 

corresponding right (left) eigenvectors |/si) and \fz) ((/si| and (/2I), a complete biorthogonal set for the spectral 
decomposition of iJr, i.e. 

Ir = \fn){h\+ l^du\U){fu\ (102) 

Ht = Ao 1/0) (/o I +J^duXu \fu){ful (103) 



ifnlh) = 1, ifulfw) = Sr{u - u'), = = 0. 



(104) 



For the H ami ltonia n Hr of the general form given in equations (101), the orthogonality and completeness relations 
(equations ( 104 ) and ( |l02|) ) can be verified by explicit calculations up to any finite order in the perturbation. For the 
Friedrichs model (a special case of the Hamiltonian given in equation (|9^)), the solutions of the eigenvalue problem 
obtained by the perturbation alghorithm coincide with the exact solutions obtained by Sudarshan C.B.Chiu, V.Gorini 
[01 . The same complex spectral decomposition for the Friedrichs model was also obtained by T.Petrosky, I.Prigogine 
and S.Tasaki |^ , with a perturbative method using a "time ordering rule" by which an small imaginary part is included 
to avoid small denominators, with a different sign according to the type of transition involved. In the perturbative 
algorithm presented in this paper, there is no need of additional rules to avoid the singularities due to resonances 
between discrete and continuous parts of the unperturbed spectrum. When states and observables are restricted to 
have suitable analytic properties, the deformation of the integrals over the continuous part of the spectrum to a curve 
in the complex lower half plane, eliminate the continuous-discrete resonances. After the expressions for the eigenvalues 
and eigenvectors are obtained up to de desired order, the complex contour of integration can be deformed back to the 
real axis. In this way, the "time ordering rule" of reference |^ can be deduced from the analytic extension properties 
of states and observables. 

If the transition probability 



Pt = |(*|exp(-zHt)|$)r 



exp(~iAot)(*|/o)(/o|$)+ / dzexp(-zzt)(*|/,)(/,|$) 



(105) 



is obtained through eigenvalues and eigenvectors computed up to n-th order, the necessary conditions for a good 
approximation of Pt are e"'^^ ^ 1 and e"^^t <^ 1 where e is the small interaction parameter. Then, even for small 
interactions, it is only for not too large times that we can obtain a good approximation of Pt. 
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An important advantage of the spectral decomposition in terms of generalized eigenvectors with complex eigenvalues 
is that it can be obtained with a well defined perturbation algorithm. The fact that this spectral decomposition is only 
feasible for states |<i>) and observables j^") with the analytic properties mentioned above is not a limitation. When 
we have to model a state or an observable from the finite amount of information coming from experimental data it 
is always possible to choose {uj\^) {{uj\'^)) with well defined analytic extensions to the lower (upper) complex half 
plane. Moreover, we are usually free to choose upper analytic extensions for the states and lower analytic extensions 
for the observables [||. For this choice, we have 



Pt = |(*|exp(-zi7i)|$} 



expHAoO(*l/o)(/n|*)+ / dz cxp{^i zt) {^\f,){f,\^) 



(106) 



where and j/^) ((/o| and {fz\) are right (left) generalized eigenvectors of being F a curve in the upper half 
of the complex plane. 

For states and obse rvable s ha ving both analytic extensions to the upper an d low er half plane near the real axis, 
both expressions (105) and (106) can be used to obtain Pt. However, equation (106) would be a bad choice for t > 0, 
as the positive imaginary part of the eigenvalues will give a well defined Pt in terms of unbounded factors (the time 
evolution would include exponentially growing terms). 

We also considered in this paper the class of observables represented by operators of the form 



\0)^Oi\l)+ dujO^\u)+ duj duj' 0^^> \ujuj') + 

Jo Jo Jo Jo 



dioO^i \ujl) + / dw'Ow |lw'), 

'0 



where the second term of the right hand side is a singular operator, and 

|1) = |1)(1|, M = \u:oj')^\u){Jl |c.l)-|^)(l|, |W) = |1)(^'I- 



(107) 



(108) 



Following the work of references |^ , Q and , we considered the states as functionals p acting on observables O to 
give the mean value of the observable in the state {{0)p = {p\0)). The states can be expressed as linear combinations 
of the functionals defined by 



(l|0) = Oi, (c^|0)=0^, (ww'IO) =0^^-, {lo1\0)^O^i, (W|0)=Ow 



(109) 



The Liouville-Von Newmann superoper ators £,= \ H, 1 and £= [H, ] can be expanded in terms of the generalized 
observables and states given in equations ( |10^ ) and (109). For the Friedrichs model we obtain 



L 



doj'{n-oj')\luj'){W\+ / duj{LU - n)\ujl){Lul\ + dw dw' {uj - uj')\ujuj'){ujuj'\, 



V^\l)+ / dw'V^^\uj'u) 



dioV^ [\ujl) - \lLo)] (1| + / dioV^ [\luj) - \loI)] {uj\ + / duj 



+ du) 



Vjl) - / dij'V^,\iJUj') 



+ I duj I duo' [V^\luj') -V^>\ujl)]{ujuj'\. 



(110) 



(111) 



1 

Assuming suitable analytic properties for p and O, the eigenvalue problem for L =L + L can be transformed into 
the eigenvalue problem for 'Lext, obtained from L replacing the integrals dio and dw' by /p and /p dz' , being 

— . . . . . 

F (F) a curve in the upper (lower) complex half plane. Starting from the right and left generalized eigenvectors of Lext, 

a well defined perturbation algorithm can be used to obtain a complete biorthogonal set of generalized eigenvectors 

of l^ext 

l^ext= j dLo\^\^^){^^\+\d\^d){'^d\+ j_du\ui\^ui){'iui\+ j du \iu'\^iu'){'^iu'\+ j_du j du' Aw | ) w | ■ 



= j + \^d){^d\ + j_du |$„i)(^'„i| + ^ du' |$i„0(*iu'| + j^c'^^ 



The generalized eigenvalues of Lext and l^ext a re sh own in Figs. 2 and 3. The time evolution of the mean value for 
any observable O of the form given in equation (107) can be obtained in terms of these generalized eigenvalues and 
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eigenvectors. For example, it is possible to obtain the time evolution of the unstable state po = If ^ 1 and 
t < (being e the interaction parameter), we obtain 

(PjI = (l|exp(iLi) ^ exp{~2TrVit){l\ + [l - exp{~2TrV^t)] {lu = n\. 

Therefore, the obtained set of generalized eigenvectors is a useful tool for a complete description of the decay process. 
The pure state |1) is given in this formalism by the state functional (1|, and the standard approximated expression 
for its survival probability is reobtained. But in addition we obtain an explicit expression for the by-products of the 
decay process. 

The generalized spectral decomposition is useful to describe the time evolution of the physical state functionals, 
because the contribution of the resonances is incorporated through the exponential dumping factors. However, it is 
important to understand that not all the generalized left eigenvectors of the Liouville-Von Newmann superoperator L 
involved in the spectral decomposition can have an independent physical meaning: Let us consider (p\ = (^'a|, where 
(^'aI is a left generalized eigenvector of L with eigenvalue A 7^ 0. Therefore we have (^'aIL = A(^'a|. Acting with 
both sides of this equation on the identity operator / we obtain 

A(*a|/) = (*a|IL/) = {-^xlHI - IH) = (*A|i^ -H)^Q. 

If A 0, ~ and the total probability condition {p\I) = 1 is not verified by (^'a|. Therefore, it can not be a 

physical state. 
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APPENDIX A: SECOND ORDER EIGENVALUES FOR THE DISCRETE SPECTRUM. 



As we showed in section II, Lext has a zero eigenvalue with infinite degeneration. This degeneration remains 
unchanged up to first order but it is partially removed in the second order approximation. 



For the right eigenvalue problem, the second order corrections A to the eigenvalues is given by equation (76) 
2TnV^ |1)(1| I) - 2mV^ \l){uj = VL\ I) =A |1)(1| !>)+ A j 
which is equivalent to the set of equations 

A (1| h = 2^iVl (1| I) - 2mV^ [u = n\ $) 



2 
A $) 



0, 



(Al) 

(A2) 
(A3) 



If \^ 0, equations (^) (|A2|) give (w| $) and A (1| $) = SttzV^J (1| $). Therefore 

\d^2mV^, \l^) = \l). 



If A= 0, equation (A3) give no condition on (w| $), so we can choose <|>) — 5{ijJ — uj). Replacing in 

_ 
(1| $) ~ 5{Vl — uj), and therefore 



2 

X~u= 0, 



I $ii) = p) + (5(<i-r!)|l). 
For the left eigenvectors, the second order correction to the eigenvalues is given by equation ( |8^ ) 



2 



2mVl (* |1)(1| - 2mVl - r!| =A (* |1)(1|+ A / dw(* \io){uj\ 



(A4) 
we obtain 

(A5) 
(A6) 



or equivalently 



20 



A(*|l) = 2^iysf(*|l), 

A i^t \uj) = -2TriV^ (^f |1) 5{lu - n). 



(A7) 
(A8) 



If At^ 0, equation (A7) gives Ad= 27rzV"Q and no condition on |1). If we choose |1) = 1, equation (A8) implies 
(I'd \lu) = -S{lu - n). Therefore 



Ad- 2TriVr 



n 1 



{■^a I = (i| - / dujS{uj - n){Lj\ = (i| -{uj^n\ 



(A9) 



If A= 0, equation (A7) impUes (ip |1) = 0, while equation (A8) gives no condition on (ip |a;). We can choose 

L-0, (AlO) 

We have obtained the eigenvalues and eigenvectors quoted in equations (|7^) and (|8^). The normalization constants 
have been chosen to satisfy the orthogonality conditions 



00 00 , 0000 

(*d|$d)-l, (^-o; = (*d I - I $d) - 0, 



(All) 



as can be easily verified using equations (A4), (|A^), ( [A9| ) and ( AlOp . It is also straighforward to verify that these 
eigenvectors form a complete set to expand the subspace generated by Pq, i.e. 



= |1)(1| + / da;\Lu){uj\ = \ $d)(*d \ + du;\ $^)(5'„ |. 



(A12) 
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